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^ . PACS. .47.54.+r- Pattern selection. 

CSl ■ PACS. .47.20.Ky- Nonlinearity. 

PACS. .42.65.Sf- Dynamics of nonlinear optical systems. 
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Abstract. - The general form of the amplitude equations for a hexagonal pattern including 
l^-^ , spatial terms is discussed. At the lowest order we obtain the phase equation for such patterns. 

f*"^ ■ The general expression of the diffusion coefficients is given and the contributions of the new 

QQ ' spatial terms are analysed in this paper. From these coefficients the phase stability regions in 

' a hexagonal pattern are determined. In the case of Benard-Marangoni instability our results 

agree qualitatively with numerical simulations performed recently. 
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Several systems out of equilibrium exhibit hexagonal patterns. Historically the cellular 
patterns reported by Benard almost a century ago were the first nonequilibrium system showing 
this planform More recently, hexagonal patterns were obtained in front solidification |2|, 
in Rayleigh-Benard convection with non-Boussinesquian effects |^ , in Faraday crispation | i | , 
' in a nonlinear Kerr medium Q , in a liquid crystal valve device , in chemical Turing patterns 
[0, in ferrofluids |^ and in vibrating granular layers Although the physical mechanism 
responsible for these patterns is different in each system they can be described within a common 
framework. A hexagonal pattern can be seen as the superposition of three systems of rolls at 
27r/3 rad, so that the resonance condition kj^ + k2 + k3 = is satisfied. The main aim of this 
paper is to discuss the form of the evolution equations of the amplitudes of those three modes, 
ie., the so-called amplitude equations, as well as the most general linear phase equation for a 
hexagonal pattern and the corresponding stability regions. 

From symmetry arguments one can deduce pO| that the normal form for the amplitude of 
the modes forming the hexagonal pattern is 

dtAi = eAi + aoMs - 7(l^2p + l^sH^i - \A,\^A, (1) 

(The overbar denotes complex conjugation. The equations for A2 and A3 are obtained by 
rotating the subindexes.) Spatial variations can be included following the Newell- Whitehead 
technique (see ref. fll|] ). Up to the third order in the amplitudes the linear spatial variations of 
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Here is the figure 1. and 2. 



Fig. 1. - Unitary vectors: parallel and ti perpendicular to the wavenumbers in a hexagonal 
pattern. 

Fig. 2. - a) Dilatations and b) distortions of a hexagonal pattern. 

each system of rolls is in the form (fi;|^ • V)^Ai, a term that must be added to eq. (|l|). Here iii 
indicates a unitary vector in the direction of the first system of rolls. Until recently, as in the 
case of a pattern of rolls, only this term was considered. Brand ||l2| discussed the possibility 
of including terms in the form (AVA). Considering the hexagonal symmetries (reflections on 
X-axis and Y-axis, and 27r/3 rad rotations) this author deduced that a term in the form 

?/3i[]43(n2- V)A2+A2(n3- V)]43] (2) 

is a real coefficient) could be also added to eq. (Q). However, as noticed by several authors 
[ p^ , p^ , in the same order another term can appear, namely 

i/52[l3(n3- V)A2+A2(n2- V)l3] (3) 

(again (32 is a real coefficient). It can easily be seen that this term also remains invariant under 
the hexagonal group transformations [p^ . 

Then, for perturbations up to the third order a generalized amplitude equation that accounts 
for spatial variations in a hexagonal pattern is 

dtAi = eAi + (nj • V)^^! + 00^2^3 + 

+ il3i\As{rv2 ■ V):42 +:42(n3 • vyA^] + i/32p3(n3 • V):42 +:42(n2 • VyA^] - 

- 7(l^2p + |A3nAi-|AipAi (4) 

The gradient terms in eqs. (|^) and (^) are a consequence of quadratic resonances and therefore 
their influence is esentially two-dimensional. (It should be noticed that, at the same order, 
terms in the form |ApVA could be included, but for the shake of simplicity we will not consider 
them in this paper). To gain some physical insight into the problem it is useful to express 
the derivatives in eq. (^) in terms of the unitary vectors of the corresponding mode, i.e, 

n2 = —563 + -^T3 in the first term and 63 = —562 — -^'?2 in the second, where stands 
for the unitary vectors perpendicular to the direction of the wavenumber of the corresponding 
system of rolls (See fig. 1). Notice that the gradient terms (eqs. (||) and (|[)) can be added to 
give 

iai[A3(n2 • V)32 +32(n3 • V)^3] + ia2[A2{% ■ -M{t2 ■ V)A2] (5) 
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with ai = /?! — 2 and 012 = -^^2- Terms in this form have been discussed by Gunaratne et al. 

[ p^Sl for chemical reactions, Kuznetsov et al. for Rayleigh-Benard convection and Bragard 

and Golovin et al. for Benard-Marangoni convection. A term with the coefficient ai 

accounts for distortions in the direction of the rolls and it therefore corresponds to dilatations 

of hexagons (it slightly changes the volume in Fourier space), while the terms with a2 account 

for distortions of the hexagonal form. In fig. 2 we represent the action of these two terms in 

Fourier space. As discussed by Kuznetsov et al. Q, ao and ai vanish when the bifurcation 

that leads to the hexagonal pattern is supercritical, while in the subcritical case the three 

coefficients a are, in general, different from zero. The term with ai will stabilize patterns with 

^ \kc\ while the term with a2 would stabilize non-equilateral hexagons. 

With the transformation in eq. (^), eq. (||) becomes 

dtAi = eAi + {hi ■ Vj^Ai + aoA2A3 + 

+ iai[A3 (n2 • V) A2 + A2 (63 ■ S/yAs] + ia2 [A2 (fa ■ \')A3 ~ A3 {h ■ V)A2] - 

- ji\A2\^ + \A3\^)Ai-\Ai\^Ai (6) 

Now let us consider solutions with a wavenumber k slightly different from kc, i.e., Ai = Aie**^' -r 
with ~ ki — kc. We are interested in homogeneous and stationary solutions of the last 
equation in the form of hexagons Ai = A2 = A3 = H ioi which the amplitude must be 

{ao + 2qai) + ^(^0 + 2gai)2 + 4(e - q^){l + 27) 

^ = WT^) 

The stability of this solution is determined considering perturbations in the form Ai — 
He^'^^'^^{l + ri + i4>i), where is the amplitude and 4>i is the phase of the perturbation. 
After introducing these perturbations in eq. (^) and linearizing one arrives at the following 
system of equations 

dtri = dlri-2qdi4)i + [ao + 2qax)H{r2+r3-ri)+H{ai + ^^{d2(t)2+d3(t)3) + 

+a2H{d3^2 + d34>3) - 2ffVi - 27i/2(r2 + rs) (8) 
dt(j)i = 2qdiri+dl(j)i~[ao + 2qax)H{4,i+<p2+(t)3)+H{ai + ^^{d2r2+d3r3) + 
2 

+ -j=a2H{d2r3 + d3r2) (9) 

where we have used the following notation 9^ = • V. We assume that the amplitudes 
and the total phase $ = (/>i + 02 + 03 are fastly decaying variables and therefore they can be 
eliminated adiabatically. As a result the dynamics is dominated by two of the phases. Instead 
of using 02 and 03 one can take a vector = [—(02 + 03), '^{^2 — 03)] and the resulting 
equation will have the most general form of a linear diffusion equation in 2D 

dt$^ Di_\7^(f+{D\\ - Dj_)\7{\7 ■ $) (10) 

This is the linear phase equation of a pattern of hexagons. The form of the coefficients in 
this equation have been chosen by analogy with that in the wave equation in an elastic solid 
[ pT| . (The velocity of the transversal waves c* corresponds here to D±, while the velocity 
for the longitudinal waves c/ is replaced by Z3||.) Using this analogy we split the phase 
into a longitudinal part 0; and a transversal 0t, that satisfy V x 0/ =0 (rhomboidal 
phase perturbations) and V • 0t = (rectangular phase perturbations). It can be proved 
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Here is the figure 3. 



Fig. 3. - Stability region of hexagonal cells with ai ~ a2 ~ 0. Notice that this region is 
symmetrical respect to the vertical axis. The phase instabilities are represented by a solid 
(£>l = 0) and a dashed line {D± = 0). 



straightforwardly that these components satisfy 

dt$i = D\^V''$u dt$t = C±vVt (11) 

A linear stability analysis of these equations shows that the system is stable to phase pertur- 
bations provided that D±^ > and Dy > 0. After tedious calculations one arrives to a general 
expression of these coefficients 

= ] - ^ + ^i^i - ^^2f (12) 



Dn 



4 2u 8u 



H Ho. Ho 
— (ai - %/3a2)^ -{ai + %/3a2) + — (3ai + ^/ia2\li) 



4 2uv 8u 

with the relationships 

u = H^{l-j) + iao + 2aiq)H >0 (14) 
V = 2H^{l + 2-f) - {ao + 2aiq)H > (15) 

The curves D± = and _D|| = determine the stability of a regular hexagonal pattern 
to rhomboid and rectangular phase perturbations respectively, while u = determines the 
region where the hexagons are unstable to amplitude perturbations. Let us mention that 
we have assumed that the homogeneous stationary solution H is positive, so ao + 2qai > 0. 
Otherwise the hexagons become unstable to a global phase change from to tt (up-hexagons to 
down- hexagons). It is interesting to examine the influence of the different parameters in phase 
stability. We first consider the particular case in which the nonlinear spatial terms are absent 
(ai = a2 = 0). This case has been considered by several authors |l^, 19, The results are 



given in fig. 3 in a {q, e) representation, where the shaded area corresponds to the stability 
region of the hexagonal pattern. Notice that due to the symmetry q —q in the diffusion 
coefficients and in u and v this figure is symmetrical with respect to the vertical axis. Near 
threshold rhomboidal phase perturbations destabilize the pattern (Z?|| = 0), but for higher 
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Here is the figure 4. 



Fig. 4. - Stability region of a hexagonal pattern when all spatial terms are included, a) Closed 
case (ofo = 1, ai = 0:2 = 0.5, 7 = 4). b) Open case (ao = 1, ai = 1, 012 = 2.5, 7 = 4). 

supercritical conditions the pattern becomes unstable by rectangular phase perturbations 
{D± = 0). These two curves intersect at the values of q that correspond to the conditions 
D± = _D|[ = 0, i.e., q — + 7)/2. (For {ai,a2) 7^ these points are quantitatively 

but not qualitatively modified.) Both curves are tangent at q = with the upper amplitude 
stability curve u = 

From eqs. (p^-(p^ one can deduce the following symmetry properties 

H{ai,q) = H{-ai,-q); u{ai, q) = u{-ai, -q); viai,q) = v{-ai,-q) (16) 
Di\{ai,a2,q) = L'||(-ai, -q;2, -g); Dj_{ai,a2,q) = Dj_{-ai,-a2, -q) 

These symmetry expressions imply that the stability curves for a particular value of (ai, 02) 
become reflected in a {q, e) representation with respect to the vertical axis under the trans- 
formation ai — > —ai, a2 — > —a2- The stability curves are qualitatively modified when 
the gradient terms are present. In general, for {ai,a2) 7^ the stability regions are no 
longer symmetrical with respect to the e-axis. An example of this situation is given in fig. 
4a. Several cases are possible. We see that for (ai,a2) > the phase instability curves 
(Z3|| = 0, D± = 0) are decentered to the right, while the minimum of amplitude instability 
curve M = is decentered to the left in the {q, e) plane. The phase instability and the amplitude 
instability curves are tangent at g = when the two conditions D± — and u ~ are met 
simultaneously. This leads to the condition ai = 1/30:2 ■ (Notice that the two phase instability 
curves are tangent at the same point). But in many cases the phase and amplitude instability 
curves intersect at two points, namely 



However, when the following condition |Q;i(ai — '\/3q!2)| > (7~ 1) is satisfied the curve D± = 
is not closed and the phase and amplitude instability curves do not intersect themselves in the 
g > quadrant. (Such a situation is represented in fig. 4b) 



ao(ai 



1 - V3a2) 



(17) 



2[ai{ai 



V3a2)) ± (7 - 1)] 
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When convection is due to surface-tension variations with the temperature (Benard-Marangoni 
(BM) convection) the condition a\ > is satisfied Two examples of the the phase stabihty 
region for {ai^a^) > are given in Fig. 3. We notice that the stabihty regions in Figs. 3a 



and 3b are in quahtative agreement with the numerical results obtained by Bestehorn |22 



for the BM problem. We shall mention that in experiments performed by Koschmieder |23 



the number of cells increases (q increases) with the supercritical heating in BM convection. 
The corresponding values are fitted quite well with the line of maximal growth rate obtained 
numerically by Bestehorn p^ . (This line remains inside the phase stability region.) This 
asymmetry could explain why in that system the transition between hexagons and rolls is not 
observed near threshold. 

From the general form of the amplitude equations we derived the phase equation for a hexag- 
onal pattern which is formally similar to the wave equation in an elastic solid. The expressions 
of the coefficients in this equation allow to determine the stability diagram. Unfortunately, 
experimental results on these phase instabilities are not yet available. We hope that the 
present results will suggest new experiments to study the wavenumber selection mechanisms 
in hexagonal patterns in different physical systems. 
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